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We investigate the role of long-lasting quantum coherence in the efficiency of energy transport 
at room temperature in Fenna-Matthews-Olson photosynthetic complexes. The excitation energy 
transfer due to the coupling of the light harvesting complex to the reaction center ( "sink" ) is analyzed 
using an effective non-Hermitian Hamiltonian. We show that, as the coupling to the reaction center 
is varied, maximal efficiency in energy transport is achieved in the vicinity of the superradiance 
transition, characterized by a segregation of the imaginary parts of the eigenvalues of the effective 
non-Hermitian Hamiltonian. Our results demonstrate that the presence of the sink (which provides 
a quasi-continuum in the energy spectrum) is the dominant effect in the energy transfer which 
takes place even in absence of a thermal bath. This approach allows one to study the effects of finite 
temperature and the effects of any coupling scheme to the reaction center. Moreover, taking into 
account a realistic electric dipole interaction, we show that the optimal distance from the reaction 
center to the Fenna-Matthews-Olson system occurs at the superradiance transition, and we show 
that this is consistent with available experimental data. 



PACS numbers: 05.50.+q, 75.10.Hk, 75.10.Pq 
I. INTRODUCTION 



The annual amount of energy humans currently use 
is delivered to Earth by the Sun in a few hours! Since 
solar energy is very dilute, it is essential to transport the 
captured energy efficiently. Most natural photosynthetic 
systems consist of antenna complexes, which capture 
photons from the Sun and transport energy to a reaction 
center (RC). There, it is transformed into chemical 
energy via charge separation. Antenna complexes are 
able to transfer excitations to RCs with an efficiency 
exceeding 95%. For a long time, it was thought that 
energy transfer in photosynthetic light-harvesting 
complexes occurs through classical processes, similar 
to random walks of the exciton to the RC. However, 
surprising evidence of coherent quantum energy transfer 
has been found recently [TJ [2] • These findings raise two 
basic questions. How can coherence be maintained in 
complex biological systems at room temperature? Why 
is quantum coherence relevant to the efficiency of energy 
transfer? 

The first question has been addressed in j3J @]- We 
consider here the second one. It is known that quantum 
coherence can speed up energy transport through a 
quantum walk, which can be faster than a classical walk 
[5]. Although the relevance of a mechanism similar 
to Dicke superradiance [5] has been also pointed out 
in [3 [7], we focus here on a different feature of the 
"superradiance transition" (ST) [10]. We show that 
ST is a dominant mechanism in an antenna complex 



described by discrete energy levels coupled to the RC, 
modeled here by a sink having a continuum energy 
spectrum similar to what has been done in [31 [51 ITT] . 
On the other hand, the effects of the thermal bath lead 
only to small corrections to the energy transport in the 
vicinity of maximal efficiency. The antenna-sink coupling 
causes the appearance of a resonance width (inverse of 
life-time) and an energy (Lamb) shift. For weak coupling 
strength, the resonance widths are roughly the same. 
However, if the coupling strength reaches a critical 
value, at which the resonance widths start to overlap, 
then a segregation of widths builds up. In this regime, 
almost the entire (summed up) decay width is allocated 
to just a few short-lived "superradiant states", while 
all other states are long-lived (and effectively decoupled 
from the environment). We call this segregation the 
"Superradiance Transition" . This effect has been studied 
using random matrix theory |12[ 113] , in nuclear physics 
[14) . for microwave billiards [15] and in paradigmatic 
models of coherent quantum transport [HI [T7] . It was 
shown in |16) that in a realistic model for quantum 
transport, maximum transmission is achieved at ST. 
In this paper, we focus on transport properties of the 
Fenna-Matthews-Olson (FMO) complex, found in green 
sulphur bacteria. This complex, one of the most studied 
in the literature [3HS1 EE] ; acts as a conductor for energy 
transport between the antenna system and the RC. The 
FMO complex is a dissipative open quantum system 
which interacts with the thermal bath provided by 
the protein environment. Here we take an effective 
non-Hermitian Hamiltonian approach [5J 1101 119] and 
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study the ST as a function of the coupling to the RC and 
the thermal bath, due to phonons. The phonon bath 
induces dephasing and dissipation, and we take both 
effects into account using two different models for the 
phonon bath. Since ST is due to quantum coherence, we 
address here two main issues: i) whether its effects can 
survive in presence of dephasing induced by the phonon 
bath at room temperature, and ii) how ST depends on 
the strength of the coupling between the FMO and the 
RC. It has been shown recently that maximal transport 
efficiency for the FMO complex is achieved near a 
critical coupling to the RC J2DJ • However, so far, the 
dependence of this critical coupling on the parameters 
of the FMO and the RC has not been determined. We 
compute this critical coupling analytically and show 
that it corresponds to the ST. We demonstrate that the 
quantum coherent effect of ST, even taking into account 
dephasing and relaxation, determines the maximal 
transport efficiency at room temperature. Indeed, 
the ST is due to coherent constructive interference 
between the various paths to the RC, thus enhancing 
the rate of energy transfer. Finally, with the aid of 
the non-Hermitian Hamiltonian approach, we consider 
a realistic coupling between the FMO complex and 
the RC, showing that the ST determines the optimal 
distance from the RC to the FMO system. 



II. THE MODEL FOR SUPERRADIANCE 
TRANSITION 

The FMO complex is a trimer, composed of identical 
subunits, each of which contains seven bacteriochloro- 
phylls (BChl) [5T]. Each subunit acts independently and 
can be modelled using a tight-binding Hamiltonian, 



(i) 



Here, |i) is the state in which the i-th site is excited 
and the others are in the ground state. Since the solar 
energy is very dilute, we limit the description to a single 
excitation in the complex, as is commonly done in the 
literature. The numerical values of Ei and Jij have been 
taken from Ref. [3J. Below we take the matrix elements 
of Hq expressed in cm~ x : 



(recombination) with an associated decay time, T\, and 
coupling to the RC with a decay time, T± r . 
As is common in quantum optics |22) . we describe this 
dissipative system with at most one excitation by states 



/ 200 -87.7 

-87.7 320 

5.5 30.8 

-5.9 8.2 

6.7 0.7 

-13.7 11.8 

V -9.9 4.3 



5.5 -5.9 6.7 -13.7 -9.9 \ 

30.8 8.2 0.7 11.8 4.3 

-53.5 -2.2 -9.6 6 

-53.5 110 -70.7 -17 -63.3 

-2.2 -70.7 270 81.1 -1.3 

-9.6 -17 81.1 420 39.7 

6 -63.3 -1.3 39.7 230 / 

The incident photon creates an electron-hole pair, called 
an exciton, which decays due to two processes: coupling 
to the electromagnetic field, i.e. emission of a photon 



M = 1>|0> ® |i> 



E 



dE b c (E)\c,E) ®\gs), (2) 



where |0) is the vacuum state of the environment and 
|c, E) ® \gs) is the state with one excitation in the envi- 
ronment and none on the sites. Here, c, is the quantum 
number labelling channels (at energies E) in the envi- 
ronments. The reduced density matrix is obtained by 
tracing over the states |0) and \c,E), 



liM+a-EN^H. (3) 



which is an 8 x 8 matrix. However, (gs\p\i) = since 
with the choice (2]), we neglect the transitions \i) — > \gs). 
Moreover, (gs\p gs) is simply the loss of probability of 
excitation of the seven sites. Therefore, we restrict our 
considerations to the 7x7 matrix (i\p\j), 1 < i, j < 7, 
which however does not have constant trace. 
In order to compute the evolution of the reduced density 
matrix, we introduce an effective non-Hermitian Hamil- 
tonian [TH1 H3] which in general can be written as, 
H cS (E) = H + A(E) - iW(E), where H is the Her- 
mitian Hamiltonian of the system decoupled from the 
environments and A(i?) and W(E) are the induced en- 
ergy shift and the dissipation, respectively. Neglecting 
the energy dependence and the energy shift we have 



H, 
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The real symmetric matrix, W, is given in terms of the 
bound state-continuum transition amplitudes, A?, from 
the discrete state i to the continuum channel c. 
The Schrodinger equation and Eq. ([2]), result in the fol- 
lowing equation for the coefficients, ay. 



iMj = 2J (Ho jk ak - iWj k a k ) 



(5) 



fc=i 



and from this the master equation easily follows, 

7 

ihp jk = [Ho , p]jk - i ^2 (Wjipik + PjiWik) . (6) 



i=i 



Under the standard assumption [2H HSJ that each site 
is coupled to an independent (local) environment, with 
associated coupling time T±, we have A] = y/h/2Ti, i — 
1, .., 7. The site i = 3 is the only one which is, in addition, 
coupled to the R C, giving rise to a decay time T lr . Then 
A| = \Jhj1T\ r (in this scheme there are 7+1 channels); 
for other channels we have A\ = 0. In Eq. @ we take 
into account: (i) the interaction between the FMO and 
the RC through the time T\ r and (ii) the characteristic 
time of exciton recombination, T±. The effects of the 



thermal bath will be considered in Sec. (IV). 
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One can verify that equation ^ can also be obtained 
by restriction to the 7x7 density matrix, from a 
Lindblad dynamics for the full 8x8 density matrix (J3|. 
In the following, we fix T\ = 1 ns, which is the exciton 
recombination time reported in the literature |24l 125) . 
and we focus on the effect of varying T± r . 



III. 



SUPERRADIANCE TRANSITION 



ST can be analyzed by studying the complex eigenval- 
ues, £ r = E r — iT r /2 of H c ff, defined in Q. As the 
coupling between the excitonic states and the RC in- 
creases, one observes a rearrangement of the widths, IV 
(the "superradiance" transition [IB] ) . We show this effect 
in Fig. [I] (left panel) , where the largest width (red dashed 
curve) and the average of the 7—1 = 6 smallest widths 
(black full curve) are plotted as functions of T\ r . For 
weak coupling to RC (large T lr ) the widths of all states 
increase as T\ r decreases. On the other hand, below a 
critical value Tf£, corresponding to ST, (vertical line), 
the average of the 6 smallest widths decreases while the 
largest width, corresponding to the superradiant state, 
increases. To examine localization of the excitation we 
use the Participation Ratio (PR)[H] of a state \tp), de- 
fined as: 

Pi?=(^|( l |^)| 4 )- 1 . 



Its value varies from 1 for fully localized to 7 for fully 
delocalized states. The right panel of Fig. [T] shows the 
PR for the state associated with the largest width (the 
one decaying most quickly). In the superradiant regime, 
Ti r < Tir, this state is fully localized on site 3, the only 
site connected to the RC. For weak coupling to the RC, 
Ti r > Tir, the PR is approximately 1.6. This small 
value, as compared to the maximal possible value of 7, 
is explained by (Anderson) localization [57] of the eigen- 
states on sites. The Anderson localization effect in the 
FMO system is due to the fact that the excitation en- 
ergies of the sites and the couplings among them are all 
different. Thus, the FMO complex can be thought of as 
a disordered system. The critical value, Tf£, at which ST 
occurs, can be estimated analytically. If all states have 
roughly the same width, at least for small coupling, then 
the superradiance condition coincides with that of over- 
lapping resonances. Such a reasoning can be applied to 
the FMO system, too. Here, eigenstates are mostly local- 
ized on the sites, and only site 3 is coupled to the RC. The 
widths are thus not uniform and most of the total width 
belongs to the eigenstate localized at site 3. Imposing 
that the half width, T3/2, is approximately equal to the 
mean level spacing D, T 3 /2 w D, and using T 3 w h/Ti r , 
we obtain the critical value at which ST occurs, 
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FIG. 1: (Color online) Left panel: average decay width, 
normalized to the mean level spacing, D as a function of the 
coupling time to the RC, T\ r . The black curve represents the 
average decay width of the 6 states with smallest width, while 
the dashed red curve shows the largest decay width. Right 
panel: PR of the eigenstate of H c n with the largest width 
as a function of T\ r . In both panels the vertical (green) line 
indicates the critical value of T\ r at which ST occurs. 



In the FMO system, the energy level spacing is 
D/hc ~ 83.5 cm' 1 , which gives T-f£ ss 0.03 ps, a value 
in very good agreement with the numerical results of 
Fig. [I] (vertical line). 

Such a value, T\ T = 0.03 ps, corresponds to a transfer 
rate, k, from site 3 to the RC of k = = 16.6 ps' 1 . 
This value is larger than the values usually mentioned in 
literature, which range from 0.25 ps" 1 to4ps _1 |25], even 
if 1 ps -1 is the most common value [5]. This discrepancy 
can be due to the simplicity of our model, even if it is 
important to notice that, to the best of our knowledge, 
the real value of the coupling time, Ti r , is not exactly 
known. In any case, for this reason, in Section VIb, we 
consider a more realistic coupling scheme between the 
FMO system and the RC. 



IV. EFFICIENCY OF ENERGY TRANSPORT 
IN PRESENCE OF A THERMAL BATH 



Interaction with the phonon environment is compli- 
cated, and it involves both dephasing and dissipation [3] . 
Since superradiance is due to quantum coherence, in Sec 



(IV A) we first focus on dephasing and the consistent 
indirect relaxation, induced by the presence of classical 
noise. On the other hand, in Sec. (IVB), we consider 



rpcr 
±lr 



h 

2D' 



(7) 



both dephasing and dissipation induced by a finite tem- 
perature bath. Needless to say, while the latter bath 
induces at equilibrium a Gibbs energy level distribution, 
the former gives rise to an equal population of all energy 
levels. 
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A. Efficiency of energy transport in presence of 



As a first step, one can study the effects of the phonon 
bath modelling the thermal bath by a classical noise. In 
this case, the dephasing effects are adequately described 
using an interaction as in Ref. [24] : 



with 



H SB =Y1 «< (*)!*><* 



(q l (t)q J {t)) = h 2 ld 6 id 6(t), 



(8) 



(9) 



where jd plays the role of the dephasing rate. This ap- 
proach corresponds to an effective infinite temperature 
that leads to equal populations of energy levels at suffi- 
ciently large times. We take into account the interaction 
^ by adding a dephasing Lindblad operator to the mas- 
ter equation ([6]), as was done in Ref. [19]. The interaction 
with noise leads to the Haken-Strobl master equation for 
the density matrix of the following form: 

*^jL = -i (H cSP - pH+ s ). ; . - ld (l - S^fHj. (10) 



The first term in the r.h.s. of Eq. ( 10 1 takes into account 



the coherent evolution and the dissipation by recombina- 
tion and trapping into the reaction center, (it is simply 
Eq. (|| rewritten in terms of H e g, defined in Eq. Q). 
The last term, which corresponds to the decay of the off- 
diagonal matrix elements, takes into account dephasing 
and indirect relaxation. For the FMO system, the ampli- 
tude of noise, hj d , is related to the temperature, T, by 
the relation, found in experiments [2], 

hj d (T) ~ 0.52 he (T/K)(cm)^ 1 , (11) 

where T is the temperature expressed in kelvin degrees, 
K. 

Transport efficiency has been measured as in Ref. [25] 
by the probability that the excitation is in the RC at the 

time tjnax: 

Wmax) = 7^- dt P33(*), 



(12) 



(13) 



Tlr Jo 

and by the average transfer time to the RC [24] 
1 f°° 

t = — / dtt pssit) Moo). 

J-lr Jo 

In our simulations, we take the initial state 
P(0) = ^(|1>(1| + |6)(6|), 

since sites 1 and 6 receive the excitation from the antenna 
system [2"4"] . 

It was numerically found in [20] that the efficiency 
reaches a maximum as a function of T lr . Here, we 
explain this as a consequence of ST, a general phe- 
nomenon in coherent quantum transport. In Fig. [2] we 
plot r\(t m ax — 5ps) (upper panel), and r (lower panel), 
as functions of T lr . The maximum efficiency of energy 
transport (maximum r\ and minimum r) is reached near 



the ST (vertical line). Note that r\{f ma x = 5ps) has a 
maximum not only in the quantum limit (^fd{T = 0.1 K), 
black dashed curve), but also considering the dephas- 
ing rate at room temperature {jd{T — 300 A'), red thick 
curve). These results show that the effects of the ST 
persist even in presence of dephasing and indirect relax- 
ation. Within the framework of the ST, the decrease in 
efficiency for large coupling to the RC, can be interpreted 
as a localization effect, see Fig. [I] (right panel). Our re- 
sults also show that dephasing can increase efficiency, 
since it counteracts quantum localization. This effect 
is known as Environment-Assisted Quantum Transport 
(ENAQT) [21 [25]. The average transfer time, see Fig. [2] 
(lower panel), has a minimum near the ST of the order 
of a few picoseconds. This time is comparable with the 
transfer times estimated in the literature. 
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FIG. 2: (Color online) Upper panel : efficiency computed at 
tmax = 5 ps, see Eq. (12 1, as a function of Ti r , for different 



dephasing rates, see Eq. (111. Lower panel: average transfer 



time, see Eq. ( 13 1, as a function of T\ r , for the same effective 



temperatures. ST has been indicated as a green vertical line. 
The initial condition is p(0) = (1/2)(|1)(1| + |6)<6|). 
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The coupling to the RC also induces a shift of the 
energy of site 3 (not only a decay width) [16 . This shift 
is assumed to be generically of the form 6 — e/Ti r , where 
e depends on the details of the coupling. We checked 
that the effect of changing e randomly, so as to produce 
up to a 50% change in the average level spacing, merely 
changes the efficiency at most by a few percent. 



B. Efficiency of energy transport in the presence of 
a finite temperature thermal bath 



In this subection, we consider the effects of energy 
transport to the RC taking into consideration the inter- 
action with a phonon bath at hnite temperature, T, as 
described in Ref. [5]. Here we consider that only site 3 
is coupled to the RC, as described above. The Lindblad- 
type master equation has the form 



dt 



= -UH cS p-pH+) +L p (p) 



v > 



(14) 



where the action of the Lindblad operator on p, L p (p), 
is described in Eq. (5) of Ref. [5]. With this choice, at 
sufficiently large time, the transition to the Gibbs distri- 
bution occurs, in absence of any other dissipative mech- 
anism, such as the presence of "sinks" . 

In Fig. [3] (upper panel), we present our results on the 
dependence of efficiency, 77, as a function of T\ r , for three 
bath temperatures. As one can see, at the ST, indicated 
as a vertical green line, the efficiency is about 0.92, and 
it weakly depends on temperature. We also mention that 
the maximal efficiency occurs near the ST. 

The value of T\ T at which one gets the maximal effi- 
ciency should not be confused with the average transfer 
time, see Eq. (13). In particular, in Fig. flower panel, 
for parameters corresponding to the ST and room tem- 
perature, Ti r 30 fs and r w 2.1 ps. 

Comparing Fig. [2] and [3j (upper panels) , one can see 
that the presence of the phonon bath significantly in- 
creases the efficiency of energy transport to the RC, al- 
most without changing the position of its maximum. We 
also would like to mention that in the presence of the 
thermal bath, the temperature effects on the efficiency 
are less significant than in presence of a classical noise, 
as in Fig. [2j Indeed dissipation helps the system to reach 
the site 3 which has the lowest energy. 

The analysis of this Section shows that the conse- 
quences of the ST are very important even in presence 
of dephasing and dissipation. For both models of ther- 
mal bath considered in this Section, the ST provides the 
maximal efficiency of energy transport. In the following 
we will consider only the model of the phonon bath pre- 
sented in Sec. (IV A), which, as shown in this Section, is 



sufficient to capture the main effects due to the phonon 
bath. 
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FIG. 3: (Color online) Upp er panel : efficiency computed 
at tmax = 5 ps, see Eq. (12 1, as a function of T\ r , for dif- 



ferent temperatures of the phonon bath. Lower panel: aver- 
age transfer time, see Eq. ( ]13[ ), as a function of T lr , for the 
same temperatures. ST has been indicated as a green verti- 
cal line. The initial condition is p(0) = (1/2)(|1)(1| + |6)(6|). 
We use as cut-off frequency of the spectral density of the 
thermal bath lj c — 150 (cm)" 1 , and a reorganization energy 
Er = 35(cm) _1 , where the latter two quantities have been 
defined in Ref. [5]. 



QUANTUM VS. CLASSICAL 



ST implies the presence of a maximum of the energy 
transport efficiency as a function of the coupling time 
to the RC, T\ r . This effect is counter-intuitive from a 
classical point of view. Indeed, the probability to escape 
(decay to the RC) for a classical particle does not de- 
crease as the escape rate (1/Ti r from site 3) is increased. 
In order to demonstrate the difference between the ef- 
fects of quantum coherence on energy transfer discussed 
above, and the corresponding classical energy transport, 
we consider a classical master equation for the population 



G 



dynamics, as in the Forster approach 
= ^2( T i,kPk - Tk,iPi) - 



>*.3 



(15) 



where Pi is the probability to be on site i, Ti t k is the 
transition matrix and the last two terms take into ac- 
count the possibility for the classical excitation to escape 
the system. The transition rates from site i to site k have 
been computed from [29] . neglecting the dependence on 
the coupling to the RC (for a classical particle, the prob- 
ability to go from any site to site 3 does not depend on 
the coupling to the RC). 
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FIG. 4: (Color online) Upper panel : quantum and clas- 
sical efficiency computed at tmax = 5 ps, as a function of 
T\ r at room temperature dephasing rate, "/d(T = 300 K), see 
Eq. (111. Lower panel: efficiency computed at t max — 5 ps, 
using Eq. (161, as a function of T\ r , at room temperature 
dephasing rate, 7d(T = 300 A"), for different coupling to RC. 
The vertical green line represents ST. As initial conditions we 
choose p(0) = (l/2)(|l)(l| + |6)(6|). 

The comparison between classical and quantum 
behavior is shown in Fig. [4] (upper panel) . The classical 
dynamics leads to a very different dependence of the 
efficiency on T\ r . Namely, the efficiency in the classical 



case does not exhibit a maximum but simply decays with 
T\ r . This shows that the ST effect is due to quantum 
coherence only. 



VI. DIFFERENT COUPLING SCHEMES 

So far, we have considered the site 3 to be the only 
one coupled to the RC. However, it is not known for 
sure which sites are connected to the RC, even though 
sites 3 and 4 are the most likely candidates, since they 
are closest to the RC [23]. As mentioned above, the non- 
Hermitian Hamiltonian formalism easily allows one to de- 
scribe different coupling schemes, which can be included 
in the effective Hamiltonian Q by properly choosing the 

coupling transition amplitude, A i , between the sites of 
the FMO complex and the RC. Note that while in the 
previous Section we indicated the channel in the RC with 
the number 3, here we label it as RC. 



A. Coupling from site 3 and 4 

Since, to the best of our knowledge, it is not exactly 
known from experimental data how the energy trans- 
port occurs, in the following we choose three different 
situations, showing that the essential features of the 
phenomenon indicated in the previous Section does not 
change too much. Specifically we consider: 

• only site 3 is coupled to RC, so that we set Ag° = 
y/h/2Ti r (as done above), 

• only site 4 is coupled to the RC, so we set A^° = 
^Jh/2T lr , 

• both sites 3 and 4 are coupled to the RC, so we set 



a: 



In a general setting the probability for the excitation to 
be in the RC at time t max = 5 ps cannot be computed 



using Eq. (12), since by merely summing that expression 
for each site connected to the RC, we neglect interference 
effects. The efficiency should be computed using 



r]'(tmax) = l-Tr(p(t max ))- 



1 



dtTr{p{t)). (16) 



Tuo 

Here, 1 — Tr(p(t max )) is the probability that the excita- 
tion leaves the system by the time t max . The last term 
in Eq. ( 16 ) is the probability that the excitation has 



been lost by recombination during this time. If there is 
just one site coupled to the RC, then Eq. (16 1 reduces 
to Eq. (pi). 



In Fig. ]4](lower panel) we show that the efficiency is 
sensitive to different coupling schemes. In particular, we 
notice that coupling through site 4 achieves a greater 
efficiency than coupling through site 3. If both sites 
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are coupled to the RC, then the efficiency is further 
improved, and the decay for small coupling times is 
smaller than that for a single coupled site. 



B. Efficiency vs position of the RC 

In the following, we consider a more realistic coupling 
scheme to RC, namely the case in which all sites of the 
FMO system have an electric dipole coupling to the same 
channel in the RC. This assumption can be justified since 
in the RC there are the same bacteriochlorophyll (BChl) 
molecules which compose the FMO system. So, it is rea- 
sonable to assume that the excitation is transferred to 
the RC by the same mechanism that operates between 
the BChl molecules in the FMO system. 

The electric dipole transition amplitude from site i to 
the RC can be written as: 
rc C r 

= 7x3 ^ ' ~ ' ^i,rc)(Prc ' Ri.RC 

(17) 

where R i RC is the distance from site i to the RC, pi is 
the dipole moment at the site i, and p RC is the dipole 
moment assigned to the RC. We take the position of 
the BChls and their dipole moments from [2U]. Here, 
we assume that the coupling strength between the sites 
and the RC is equal to that between the sites, so that: 
C\p\ 2 = 134000 cm" 1 (I) 3 , as in Ref. [2U]. 

In order to determine the coupling amplitude from site 
i to the continuum of states in the RC we evaluate the 
transfer rate, Hi from the Fermi-golden rule [5]: 



2 ^P RC \ V i 



(18) 



where p RC represents the density of states in the RC. It is 
interesting to observe that an expression for the transfer 
rate similar to ( 18 ) can be also obtained without per- 
turbation theory, see for instance Ref. [3D], where the 
continuum was modeled as a semi-infinite lead. We can 
now determine the non-Hermitian Hamiltonian, Eq. Q, 
setting: 

, RC i RC 

A = V 2 ^Prc Vi . 

Note that now the coupling between the FMO complex 
and the RC depends on the position of the RC. In or- 
der to determine how the efficiency of energy transfer 
depends on the position of the RC w.r.t. the FMO com- 
plex, we assume p RC = pz and place the RC in the same 
y and z position of site 3, as given in [20]. The transport 
efficiency has been studied by varying the distance, d, 
from the site 3 along the x direction. In order to com- 
pute the transition amplitude, A i , we need the density 
of states of the RC which, to the best of our knowledge 
is not known experimentally. For this reason we consider 
different densities of states, respectively larger, equal or 
smaller than the density of states of the FMO system, 



p FMO ~ 1 /D, D being the mean level spacing of the FMO 
complex discussed in Sec. (III). We show in Fig. ^ how 




5: (Color online) Efficiency, Eq. (161, computed at 



FIG. 

tmax ~ 5 ps vs the distance d (nm) of the RC from site 3 in the 
FMO complex for different ratios x = Prc I Pfmo- H ere i 811 
electric dipole coupling between the RC and the FMO system, 
see Eq. (171, has been considered. The vertical dashed lines, 
obtained from Eq. (19 1, represent the critical distances at 
which superradiance transition occurs. As initial conditions 
we choose p(0) = (l/2)(|l)(l| + |6)(6|). Data in this figure 
refer to a room temperature dephasing rate, 7<j(T = 300 A'), 
see Eq. (flil. 



the efficiency, computed with Eq. (|16| using Eq. (10 1, 



varies as a function of the distance from the RC to site 
3 of the FMO system. In Fig. Q we consider a room 
temperature dephasing rate 7d(T = 300 AT), for different 
ratios, x — Pro/Pfmo = 0.1, 1, 10 , respectively, from 
the upper to the lower panel. As one can see, the opti- 
mal distance, which is the distance that maximizes the 
efficiency, slowly depends on the density of states in the 
RC. 

We can use the superradiant criterium obtained in 
Eq. |7]) to get an analytical expression for the optimal 
distance of the RC from the FMO complex. Since site 3 
is the closest to the RC, we can use Eq. ^ and we find 
that the ST occurs for Hk^ = D. Finally, from Eq. (18) 
we have: 



(2ttB 2 )Hp rcPfmo ) 



where, 



B = C 



T 



R., 



(19) 



) , (20) 



see Eq. (17). Eq. (19) gives us the distance at which 
superradiance transition occurs. The critical distance 



obtained from Eq. ( 19 1 is shown in Fig. ([5| as dashed 
vertical lines. As one can see, the estimate is very good. 
Note that changing the density of states in the RC from 
p FMO /10 to 10 2 p FMO only changes the optimal distance 
from 1 nm to 3 nm. These distances are consistent with 
available structural data for the RC-FMO complex, see 
for instance [3T]. The result is remarkable since it shows 
that the superradiant criterium suffices to determine the 
optimal distance from the RC to the FMO complex for a 
wide range of values for the density of states of the RC. 



VII. CONCLUSION 

We have analyzed energy transport in the FMO system 
with the aid of a non-Hermitian Hamiltonian approach. 
This allows us to take into account the effect of the cou- 
pling of the FMO system to the reaction center in a con- 
sistent way, not merely phenomenologically, as is usually 
done in the literature. We have shown that by increas- 
ing the strength of the coupling to the reaction center, 
a superradiance transition occurs. This transition oc- 
curs at approximately the same value of the coupling for 
which energy transport efficiency is maximal. Indeed, the 
superradiance transition is due to coherent constructive 
interference between the paths to the RC, and this effect 
enhances the rate of energy transfer. Since the ST effect 
is due to quantum coherence, one might expect that any 
consequences of ST would disappear in the presence of 
dephasing and relaxation provided by the thermal bath. 
On the contrary, we have shown that the effect of su- 
perradiance survives in the presence of the thermal bath, 
and the maximal efficiency only depends weakly on tem- 



perature. We have also estimated the ST critical value 
analytically. For coupling strengths of the FMO system 
to the RC near the critical one, where the superradi- 
ance transition takes place, we obtained average energy 
transfer times comparable to experimental values (a few 
picoseconds). Finally, we took into account a realistic 
coupling scheme between the FMO system and the RC, 
and derived from the superradiance condition the analyt- 



ical expression, Eq. (19), for the optimal distance from 
the RC to the FMO complex. This analytical expres- 
sion depends on the density of states in the RC. Within 
a wide range of density of states the optimal distance 
which we obtained analytically is approximately a few 
nanometers and is consistent with available structural 



data on RC. Note also that Eq. ( 19 ) is valid for a generic 
Donor- Acceptor complex. 

Our analysis shows that the superradiance mechanism 
might play an important role in explaining the efficiency 
of quantum transport in photosynthetic light-harvesting 
systems and in engineering artificial light harvesting sys- 
tems. 
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